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Abstract 

We study the one- and two-holon excitations of the SU(3) Kuramoto-Yokoyama model on the 
level of explicit wave functions, and generalize the calculations to the case of SU(n). We obtain 
the exact energies and the single holon momenta, which we find fractionally spaced according to 
fractional statistics with statistical parameter g = 1/n. 
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I. INTRODUCTION 

Since its discovery in 1988 by Haldane [jj and Shastry Q], the Haldane-Shastry model 
(HSM) has amply contributed to our understanding of fractional quantization in one- 
dimensional spin chains. The model provides a framework to formulate and analyze spinons, 
the elementary excitations of one-dimensional spin chains, at the level of explicit wave func- 
tions Q]. In particular, it was realized through this model that spinons in SU(2) spin 
chains obey half- Fermi statistics Q]. Kawakami y] subsequently generalized the HSM from 
SU(2) spins to SU(n), a model in which the spinon excitations obey fractional statistics with 
statistical parameter (1 - l/n) QBEHlQ, 11, 12, 13 1 - 



The HSM was also generalized by Kuramoto and Yokoyama to allow for mobile holes. 
The Kuramoto- Yokoyama Model (KYM) hence contains spin and charge degrees of freedom, 
described by spinon and holon excitations [15], which carry spin | but no charge and charge 
+ 1 but no spin, respectively. While explicit wave functions for one- holon states in the 
SU(2) KYM were known for many years jlfjj], the construction of the exact two-holon states 
was achieved only recently In particular, the single-holon momenta in these states 

were found to be shifted by a fraction of the units 2ir/N appropriate for a chain with TV 
sites, periodic boundary conditions (PBCs), and a lattice constant set to unity. This result 
was interpreted as a manifestation of half-Fermi and hence fractional statistics among the 
holon excitations [3], thus confirming a conclusion reached by Ha and Haldane ^| using 
the asymptotic Bethe ansatz, by Kuramoto and Kato [3, 0] from thermodynamics, and 
by Arikawa, Yamamoto, Saiga, and Kuramoto Q Q from the electron addition spectral 
function of the model. Like the HSM, the KYM can be generalized to spin symmetry 
SU(n) HQ- 

In this article, we analyze the one-holon and two-holon excitations of the SU(n) KYM 
on the level of explicit wave functions. The article is organized as follows. In Section |Hj we 
investigate the case of SU(3). We first present the basic properties of the model including 
the ground state and the coloron excitations in the absence of holes, where the SU(3) KYM 

n 

reduces to the SU(3) HSM studied previously in a similar framework [L3J. We then construct 
the explicit one-holon and two-holon wave functions and derive the exact energies and single- 
holon momenta. In Section ITTT1 we generalize the results to SU(n). In particular, we review 
the basic properties of the ground state and the SU(n) spinon excitations before we derive 
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the one-holon and two-holon wave functions including their energies and momenta. In 
Section ITVl we interpret our results in terms of free holons obeying fractional statistics with 
statistical parameter g — 1/n. 



II. SU(3) KURAMOTO YOKOYAMA MODEL 
A. Hamiltonian 

The SU(3) Kuramoto-Yokoyama model (KYM) 

[6] is most conveniently formulated by 
embedding the one-dimensional chain with periodic boundary conditions into the complex 
plane by mapping it onto the unit circle with the sites located at the complex positions 
r] a = exp(i^a), where N denotes the number of sites and a = l,...,N. For the SU(3) 
case, the sites can be either singly occupied by a fermion with SU(3) spin or empty. The 
Hamiltonian is given by 

7T 2 N P 

#SU(3) = — i — fo) (!) 

N a^i yn^-m 

where P a p exchanges the configurations on the sites rj a and r/p including a minus sign if 
both are fermionic. Rewriting (JTJ in terms of spin and fermion creation and annihilation 
operators yields 

2 N 

L. yj _l 

(2) 

where we label the SU(3) spin or color index a by the colors blue (b), red (r), and green (g). 
The Gutzwiller projector enforces at most single occupancy on all sites, and is explicitly 
given by 

N 

Pq = I I (1 ri a i,n ar ^ab^«g n ar n ag ~\~ ^ n ab n ar n ag) > (3) 



27r 2 1 „ 1 \ - / t t \ n a n 1 

H su(3) = ^2^ \ Vry _ Va \2 P v "2 1^ {<^p a + c 1 p a c aa j+J a -Jp-— + n a -~ 



o-=b,r„t 



Pi 



G • 



a=l 



where n a = c^c^+c^c^+c^c is the charge occupation operator at site rj a . Furthermore, 
we have introduced J a — \ J2ar c ao-^o-TC aT , the eight- dimensional SU(3) spin vector, where 
A denotes the vector consisting of the eight Gell-Mann matrices (see App. [X}, and a and 
t are again SU(3) color indices. For all practical purposes, it is convenient to express the 
SU(3) spin operators in terms of colorflip operators e a J = <^ aa c ar . The Hamiltonian (J2J) then 
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becomes 

9^r 2 JL 1 f 1 / \ 1 n n„ 1 

Pg, 
(4) 



#SU(3) - "Tp 



2 ( C «0- C /3<T + C /3(7 C C«T J + 2 e ° 



2 " 2 

cr,r 



where the color double sum includes terms with a = r. 

The KYM is supersymmetric, i.e., the Hamiltonian ([!} commutes with the operators 
Q a aa a a fe5 where a aa denotes the annihilation operator of a particle of species a (a 
runs over color indices as well as empty site) at site rj a . The traceless parts of the operators 
jab g enera te the Lie superalgebra su(l|3), which includes in particular the total spin opera- 
tors J = J2a=i J a- m addition, the KYM possesses a super- Yangian symmetry which 
causes its amenability to rather explicit solution. 



B. Vacuum state 



We first review the state containing no excitations, i.e., neither colorons nor holons. This 
vacuum state is the ground state at one third filling, where the SU(3) KYM reduces to 
the SU(3) HSM. The vacuum state for N = 3M (M integer) is constructed by Gutzwiller 
projection of a filled band (or Slater determinant (SD) state) containing a total of N SU(3) 
particles obeying Fermi statistics 

|*o> = Pg II 4i4 1°) = P g |*sd>- (5) 

l?l<9F 

As h&sb) is an SU(3) singlet by construction and Pq commutes with SU(3) rotations, |^o) 
is an SU(3) singlet as well. 

If one interprets the state |0 g ) = Yla=i c <L g I ) as a reference state and the colorflip 
operators e bg and e rg as "particle creation operators" , the state (jSJ) can be rewritten as j^Ql 

|* ) = *«»[*; v,k] e h z f • • • • • • |0 g ), (6) 

where the sum extends over all possible ways to distribute the positions of the blue particles 
zi, . . . , z Ml and red particles w±, . . . , wm 2 ° ver the N sites. The vacuum state wave function 
is given by 

Mi M 2 Mi M 2 Mi M 2 

^o[zu W h ) = - Zjf \\{w k - W{f Y[ Y\_( z i ~ w k) \\ Zi W k (7) 

i<j k<l i=l k=l i=l k=l 
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with Mi = M2 = M, its energy is 

36 V N J V ; 

The total momentum, as defined through e tp = ^o[r]iZi, r)iWk]/^?a[zi, Wk] with 771 = exp(z^), 
is p = regardless of M. For further purposes, it is important to note that the wave function 
(J7J) can be equally expressed by any two sets of color variables, as it is shown in App. O 



C. Coloron excitations 



Let N = 3M - 1, Mi = (N - 2)/3, M 2 = (N + l)/3. A localized coloron at site "?? 7 " 
is constructed by annihilation of a particle with color a from a Slater determinant state of 
iV + 1 fermions before Gutzwiller projection jllll ]: 

|^> = ^ G c 7(T K +1 >, (9) 

where a denotes the complementary color of the coloron. The annihilation of the fermion 
causes an inhomogeneity in the SU(3) spin and charge degree of freedom. The projection, 
however, smoothes out the inhomogeneity in the charge degrees of freedom, the coloron thus 
possesses color, but no charge. The wave function of a localized, e.g. anti-blue or yellow, 
coloron is given by 

Ah 

%{zi, w k ) = JJ(r/ 7 - 3i)tfofe; Wk], (10) 
i=i 

with \l/o as stated in ((Zj). Fourier transformation yields the momentum eigenstates 

1 - 

TO = ]v5>r) B |*$>' (11) 

7=1 

which identically vanish unless < n < M\. In particular, this implies that the localized 
one-coloron states (jHJ) form an overcomplete set. It is hence not possible to interpret the 
"coordinate" 7/ 7 literally as the position of the coloron. The momentum of (jllj) is 

Pn = y " ^ [n + 3 J , < n < Mi. (12) 

The momentum eigenstates (fTTj) are found to be exact energy eigenstates of the Hamilto- 
nian with energies 

K = Eo + l^- 2 +e c {pl), (13) 



where the one-coloron dispersion is given by 

e c (p) = j(y-(p-7r) 2 ). (14) 

Colorons obey fractional statistics, the statistical parameter between color-polarized colorons 
is given by g = 2/3. 

D. One-Holon excitations 

1. One-holon wave functions 

If we dope holes into the SU(3) spin chain, this will cause the existence of ho Ions, the 
elementary charge excitations of the system. In this section, we will construct the wave 
functions of the one-holon states and prove by explicit calculation that these states are 
eigenstates of the Hamilitonian (0). For this consider a chain with N = 3M + 1 sites. A 
localized holon at lattice site ^ is constructed as 

l^°>= c ^o4 CT |^ D 1 >, (15) 

where the color index a can be chosen arbitrarily. Compared to the coloron, we eliminate 
the inhomogeneity in color while creating an inhomogeneity in the charge distribution after 
Gutzwiller projection. Thus the holon has no color but charge e > (as the charge at site 
T)£ is removed). Note that the holon is constructed as apparently being strictly localized at 
the coordinate £, as states (|15j) on neighboring coordinates are orthogonal. In total, there 
are N independent states of the form (|15|) . 

Momentum eigenstates are constructed from by Fourier transformation. We will 
show below that only (N + 5)/3 of them are energy eigenstates, and restrict ourselves to 
this subset in the following. In order to describe these states by their wave functions, we 
take 1 g ) = ria=i c ag |0) as reference state and and write the one-holon states as 

= £ *m [*; h] c hg e h z f . . . . . . e^ M2 |0 g ) , (16) 

{Zi;w k ;h} 

where the sum extends over all possible ways to distribute the blue coordinates Zi, the 
red coordinates w^, and the holon coordinate h over the N sites subject to the restriction 
h 7^ Zi,Wk- The one-holon wave function is given by 

Mi M 2 

*%[zii w k ; h] = h m ]J(h - zt) ]J(h - w k )^ [z t ; w k }. (17) 

i=l k=l 
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To increase readability of the following calculations, we will keep the distinction between 
Mi and M2, although we will always set M\, M2 and M3, i.e., the numbers of blue, red, 
and green particles, to be equal to M at the end. In order for (|17j) to represent energy 
eigenstates, the integer m has to be restricted to 

< m < M + 1 = N + 2 . (18) 

3 

For other values of m, the states are not eigenstates of the Hamiltonian ([Tjl. although 

they do not vanish identically (as the l^^'s do). Consequently, we are allowed to refer to 
the states (JUJ) with (JT2J) as "holons" only ifO<m<M + l. 

This also implies that the states ()15|) do not really constitute "holons" localized in position 
space, but only basis states which can be used to construct holons if the momentum is chosen 
adequately. Since the states (fTS")) are orthogonal for different lattice positions £, there are 
N = 3M + 1 orthogonal position basis states j^ )- These states cannot strictly be holons, 
but rather constitute incoherent superpositions of holons and other states. It is hence not 
possible to localize a holon onto a single lattice site. The best we can do is to take a Fourier 
transform of the exact eigenstates hi^) f° r < m < M + 1 back into position space. The 
resulting "localized" holon states will be true holons but will not be localized strictly onto 
lattice sites. 

The momentum of ()16l) is 

The one-holon energies are derived below to be 

Et° = E -l^- 2 +e ho (p%), (20) 
where the one-holon dispersion is given by 

6 (p) = -4 ( y - (p - *) )> y- p ~Y- (21) 

In the following subsection we will prove that the states (fTT)|) are energy eigenstates of the 
Hamiltonian if (and only if) the momentum quantum number m is restricted to |T%|) . 

2. Derivation of the one-holon energies 

To evaluate the action of i?su(3) 011 |^m)' we ^ rs ^ re pl ace e T e p g h a — e^ b — e")(l — 

hp — e h p — e^ r ), where h a denotes the hole occupation operator h a = 1 — n a , and rewrite the 
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Hamiltonian (HI) as 



Hi 



SU(3) 



) 



2?F 1 f P b g gb , rg gr , br rb 

2tt 2 * 1 

] V^ 1 fpbbpbb , rr rr , bb rr\ 

^ N 2l^\ |2 l e « e /J + e «^+ e Q e /3j 



2vr 2 



AT2 



^ 1 2 2 ^ 1 / 1 



.2 ^ 



27T 2 \ 1 1 , + + . 1 . + j. 



1 , ■]- ■]- , 



(22) 



In the following we evaluate each term of (j22j) separately. 

The first term [e^ s el [zf, Wp.', h], which vanishes unless one of the z^s is equal to 
r] a , yields through Taylor expansion (the derivative operators are understood to act on the 
analytic extension of the wave function) 

Mi N tW , , 

rjp v™[...,Zi- 1 ,T}p,z i+1 ,...;w k ;h\ 



pbg p gb 

a 13 \]> ho 



2 m 



i=i |j '^ ,| 



? //3 



i=l ^=0 



94 



^{N 2 + 8M 2 - 6Mi(iV + 1) + 3) ^ 



ho 
fit 



M 2 

N — 3 \ \ - z 



Mi 



ho 



i=l fc=l 
A/i A/ 2 



3? 



ho 



Afi Af 2 



i= 
Mi 

+£ 



2 ^t^ (z, - W k )(Zi - w t ) 



ho 



2z? 



(zi - Zj)(zi - h) 



AT q M l 

_ V - ^ ho 

2 t^Zi-h m 

i=i * 



(23) 
(24) 
(25) 

(26) 

(27) 

(28) 
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Mi M 2 2 

+gg (Zi - rot V^ (») 

where we have used deg z .^ r ^[^; w^; /i] = A" — 1 and defined Ae = — ^2 a =l VaiVa — l) e ~ 2 . 
Evaluation of the latter yields A = (N - l)(N - 5)/12, A x = -(N - 3)/2, A 2 = 1, and 
Ae = for 2 < £ < iV — 1 (see App. \D§. Furthermore, we have used 

+ 7 T7 7 + 7 77 T = 1 > x,y,zeC. (30) 



(a;-y)(ar-z) (y - x)(y - z) (z-x)(z-y) 
The second term [e^fe?^^] w^; h] can be treated in the same way, yielding together 
with the first term in (|23|) 

AT „ A/i M 2 AT Q Mi M 2 , r Q 

A — 3 \ - x -v Za A — 3 \ - x w k A — 3 , , , , 

> > — + > > — = M X M 2 . 

2 ^ ^ Zi-w k 2 ^ ^ Zi - w k 2 
i=i k=i 1 K i=i it=i * K 

One part of (|26j) and the term corresponding to (}2~Tj) can be simplified with (j3*Uj) to 

Mi M 2 / 2 i 2 \ i 

7 3 T + -7 S r J =-M 1 (M 1 -l)M 2 , 

~r/ ^ V I 2 * _ *j) - 2 (^ - w fc ) (z,- - w fc ) J 2 

as well as similar expressions for z\ <-» w k . 

The third term [e^ r e^ b \I/^] [zi, w k ; h] leads to 

P br„rb 

C " C /3 lT ,ho 



Mi M 2 Mi / \ M 2 / x 



Mi M 2 

yy ZiWk ^ (31) 

«=i fc=i v j 

Mi Af 2 Mi M 2 

Mi M 2 Mj-1 , v 2 

+EEEiEr i ^ — t*S (33) 

1=1 fc=l /i=2 ^ { aj } V ai * y V ^ V 

Mi M 2 Mi-1M 2 -1 (]\v 

+ E E n\ u \ 

i=l k=l /x=l y=l P " 



{tax - Zi)--' {z a - Zi)(w bl ~W k )--- (w bv - W k ) 
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where {aj} {{bi}) is a set of integers between 1 and Mi (M2). The summations run over 
all possible ways to distribute the z a . (u>&J over the blue (red) coordinates, where Zi (u>&) is 
excluded. The two terms and vanish due to 13] 

Theorem 1 Let M > 3, z 6 C, and z 1; . . . , 2;^ G C distinct. Then, 



M 



\M-3 



E m _ = 0- (36) 

i=l Ylj^i( Z j Z i) 

The last term (J35j) can be simplified using a theorem due to Ha and Haldane (2^]: 

Theorem 2 Let {aj} be a set of distinct integers between 1 and M\, and {bi} a set of 
distinct integers between 1 and M 2 . Then, 

EEE E E - ' " A- 1 - ' " ' 



i=lfc=1 M=1 , =1 { } a*m 



min(Mi,M 2 ) 



K=l 



Furthermore, the two terms in line (JH2J), together with the remainder of (|2T)j) and the 
corresponding expression from the second term of the Hamiltonian, can be simplified to 
M 1 M 2 (M 1 + M 2 - 2)^/2. 

The 2nd and 3rd line of (J22j) yield 



N p bb bb 1 rr rr _i_ bb rr _ bb _ it , _ 1 

V-> e * e (3 + e a e (3 + e a e (3 e a e * + ™a 2 lT ,hor . . , i 

lla //31 

1 Mi 2 M 2 2 

= -(Mi (Mi - 1) + M 2 (M 2 - " E 77371a*S - E fti; !* ^ S? 

M i M 2 AT 2 / AT \ 

- E E (^«5 - ^ + * - t + • (^) 

by which the remainder of ([25)1 . its counterpart from the second term, and (|31|) are cancelled. 
The 4th line of © yields 



E (e °V eg)(1 , 2 n> ' ) ^ *] = (E + E *£• (as) 

We will now evaluate the charge kinetic terms, which include the technical improvements 
compared to previous calculations. We will use a Taylor expansion as in (|2*3Jl . For the 
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treatment of the charge kinetic terms it is crucial that the fermionic creation and annihilation 
operators appearing in the expansion match with the variables of the analytically extended 
wave function. We wish to stress that the one-holon wave function (J 17)) can, as the ground 
state wave function, be equally expressed by an arbitrary pair of sets of color variables. In 
(|22|) we have thus written the charge kinetic terms in a symmetricized way. For the first 
term we get 



N 1 t i t 

J- C ab C (3b ' C nr C 



E 



a* Pt ^ho 



\Va ~ Vp\ 



Zi] w k ; h] 



' N 

E 

&*P 



C av C f3v 
Va-VP 



_\]) ho 

2 rn 



[v u . . .,v Ml+M2 ; h] 



N 

E 

P+h 



Vp ^^[v 1: ...,v Ml +M 2 ;r]p\ 



h-Vp\' 



Vp 



Mj+M 2 N 



E 



VpiVp - hf & ( tf£fa, . . . , v Ml+M2 ; Vl3 ] ' 



^ ^ £\\h 



Vp\ 2 



M1+M2 ry m , m+ e n£ 



E 

1=0 

N 2 



12 

iV 2 -l 



£\ dh e 

1 m(m ■ 
- + — ^— r 



r lp 



V0=h 



% 2 ; h\ 



N) 



h m - 



1 

+ 2 



12 

'Mi 

E 



+ 



mm 



h m 

N -1 
2 

iV- 1 



1 «9/i + 2 dh 2 



V%[zi',w k ]h] 



h< 



h 2 



Mi M 2 



(h - Zi)(h - Zj) 



+ EE 



— m 



2/i 2 



'Mi 

E 

i=i 



M 2 

E 



M 2 

+ E 

fc=i 



h 



ho 



/r 



^ (h - Zi)(h - w k ) ^ (h - w k )(h - wi) 



*m,(39) 



L«^i v ' " /v " v ' i=i fc=i v ' " /v " """" fc^z 
where the v^s denote the union of the blue and the red coordinates, and we have introduced 
B ? = - J2a=i Va +1 (v* ~ l) e ~ 2 - Evaluation of the latter yields B™ = (N 2 - 1) /12 + m{m - 



N)/2, By = m-(N - l)/2, Bf = 1, and B e = for 3 < £ and < m < (JV + 2)/3 (see 
App. [DJ). The restriction (|18|) of the allowed momentum values follows from the 5-series in 
(|39|) . since 5^ 7^ for 3 < £ and (N + 2)/3 < m, in which case the calculations above are 
not feasible anymore. 

For the evaluation of the remaining two charge kinetic terms we reexpress the wave 
function by the other pairs of sets of color variables (see App. ICj) . We then proceed as 
in ([39)1 . where we replace the green variables by the blue and red ones using the identities 
of App. Doing so, we finally arrive at 



' N t 1 c J 

El C ab C /3b ^ C Qg C /3g , T ,ho 
2 I- -12 m 

!X+P 



Zi, w k ; h] 
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12 



A/2 



+- 



M2 



M 2 



/ — ' h — w k 



k=l 



■^{h- w k ) 



fe=i 

M 2 

E 



h-w k 



ho 



2 ^ ( h - w k)(h - wi) 



*m,(40) 



as well as 



El CcM-Cft. + C ag C ffg ho 
9 1 



2 |?7a-^| 



AT2_i m(m-iV)^ ho 



12 



1 

+ 2 



Mi 



C*i — C*2 — 2C*i ^ 



h 



N — 1 
2 

Mi 



Mi 



i=i 



ft - Zi 



ho 



Mi 



i=l 



h - Zi 



-^{h-Zif ^ (h - Zi)(h - Zj) 



*£■ (41) 



In (gUJ) and flU) we have denned the constants C x = Ea=i V(l ~ = (N - l)/2 and 
C 2 = J2a=i V(l " ^«) 2 = -(iV 2 - 6iV + 5)/12 (see App. 0. 

Now there occur several simplifications. The first term in (|2*H|) together with the appro- 
priate terms in (jHO]) and (}4*T|) yield 



2z? 



Mi 



/> 2 



Mi 

(*• - - h) ' j^,{h- Zi)(h - zj) 



Mi (Mi - 1) 



(42) 



and the similar expression for Z{ <-> leads to M2(M2 — 1). Furthermore, (j29jl . its coun- 
terpart from the second term, and the appropriate term in result in 



Mi M 2 

EE 



Mi M2 

+EE 



Mi M 2 



1 ^ (zi-w k )( Zi -h) ^^[{w k - Zi )(w k -h) j^fr[(h-z l )(h-w k ) 



+ EE 



MiM 2 . 
(43) 

Finally, the remainder of 1)280 . the first term of ()38)h and the remaining off-diagonal terms 
of (jUJ) yield 



A* — 3 2j 1 

2~ ^ Zi -h + ^ \zi-h\< 
i=i 1 i=i 1 * 1 



Mi 



N-l^K h 



Mi 



/i 2 



2 Z — ' ft — Zi 



•[ (h — Zi 



N - 3 
-Mi—— , (44) 



and the similar expression for z^ ^ w k . 
Summing up all terms, we obtain 



#SU(3) |*m) = £ m |*m) 



(45) 



12 



with 



E r , 



where we have set Mi 
into the form (|20p. 



" w 

M 2 



72 24 18 



3 / iV + 2 

—TO TO 

2 V 3 



(46) 



M = (iV — l)/3. Using (fT9*j) . can now be easily brought 



E. Two-Holon excitations 

1. Momentum eigenstates 

We will now investigate the two-holon eigenstates. For this, let the number of sites be 
given by N = 3M + 2. The state with two localized holons is constructed as 

l^6> = c a^Gi.iJ^D 2 >- (47) 

Similar to the one-holon case, these localized states (|47jl do not really represent "holons" 
localized in position space, and we can refer to true physical holons only in momentum 
space. The two-holon momentum eigenstates will be most easily described by their wave 
functions 

w k ; h, h 2 ] =(hi - h 2 ){h™h n 2 + hlh™) 

Mi M 2 (48) 

i=l fe=l 

where h 12 denote the holon coordinates and the integers to and n are restricted to 

N +1 

0<n<m<M + 1 = . (49) 

3 

This restriction will be derived below. 

The two-holon state represented by (}4~%j) is 

= £ ^ ^ <ww£ • • • « • • • e -M 2 log) > (so) 

{2 i ;w ft ;/ii,ft,2} 

where the sum contains the restriction h\ 2 ^ w^. The total momentum of the states (J5(J|) 
is found to be 

hn 4vr 2tt / l\ 

Pm« = Y+ jy ( m + 7l -3j mod27r - ( 51 ) 

In the following two subsections we construct the two-holon energy eigenstates starting 
from (g7j). T 
SU(2) KYM 



re used strategy is similar to the construction of the two-holon states in the 

3- 
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2. Action of Hgurg) on the momentum eigenstates 



In order to derive the action of the Hamiltonian on the momentum eigenstates (|50|). we 
first define the auxiliary wave functions 

Mi M 2 

(p mn [zi;w k ;hi,h 2 ] = h^h^Y^ih - Zi)(h 2 ~ z^Y^ih - w k )(h 2 - w k )^ [zi-,w k } 

i=l fc=l 

= iJ hlh2 ^ [zi;w k }, (52) 
which can be used to express the wave functions (J4~HJ) as 

^mn = ¥m+l,n + <Pn+l,m ~ V^m.n+l — <pn,m+l- (53) 

In agreement to the one-holon case, we use (|22j) for the Hamiltonian and concentrate on 
the terms which differ from the ones above. The first term [e^ g ejgV mn ] [zf, w k ;hx, h 2 ] yields 



V e h Je eh 



Mi N-l 

E 



a °/3 

2^rnn 



[zi;w k ;hi,h 2 ] = EE 



A-eZi a (p r< 



t\ dzf Zi 

i=l £=0 1 1 



^-(N 2 + 8M 2 - 6M 1 (N + 1) + 3) <p. 



N o Mi M 2 Mi 2 

4=1 fc=l l^J v J/ 

Mi M 2 2 A/l ^2 2 

+ 2 E E 7 § n + ? E E 7 f? (55) 

^fl d 2 ^ 2z 2 d N-3 d \ , 

1=1 \ 1 j^i J / 

Mi M 2 2 a 

Zj — W k OZi 

i=l k=l 

Now, the lines (p>fj) and (jB3j) can be treated as in the one-holon calculation, whereas the 
lines (}56|) and (|H7j) explicitly yield using (J30|) : 

Ml f h 2 h 2 \ 

^ t 1 " (S - - " S - ^l)(^2 - *i) ^ 



1=1 



7JWI 



fri ~ z i J 2 V h 2 — Zi 



mil 
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Mi M 2 , 

+EE> 

i=i k=i v 

Mi M 2 , 

+EE' 

i=i fc=i v 



(w k - Zi)(w k - hi) {h x - Zi)(hx - w k ) 



vol 



h\ 



(58) 



(w fe - Zi) (w k - h 2 ) (h 2 - Zi)(h 2 - w k ) 

The term [e^ g e|Vmn] [zf, w k ;hi, h 2 ] leads to the analog result with zi and w k interchanged. 
Furthermore, the terms [ejj'e*^] [z{, w k ; hi, h 2 ] as well as the 2nd and 3rd line of are 
unchanged as compared to the one-holon case. The 4th line of (|22j) yields 

Mi M 2 



N bb I rr / f m l ] JW 2 | 

E^dsd-^)^^ Ep-^ + E^^ + ^^J I - <■'<<» 

VI. i=l 1 1 k=l ' 

where {hi «-> /z 2 } denotes the reappearance of the preceeding terms in curly brackets with 
hi and h 2 interchanged. 

For the charge kinetic terms we obtain in analogy to the one-holon states 



N 



9 \Z I~75 ^mr, 



2 ka-^ 2 



N 



E l. ,^12 ^ 



,a=h\,h 2 /3^« 



- ^7/3 1 



[i;*; /2 2 ] 



Mi+M 2 nmh m+£ o£ / 

I\ dhj \~W 



M1+M2 jjnin+l ^ / 

EDg n 2 O ( (fmn 



dh e 2 V K 



N 2 ~l + m{m - N) + n(n-N) ^ ^ 
N - 1 



2 

N 2 -l 
N-l 



1 d 



1 



d 2 



1 d 2 

\min+2 



— n I h m h n+l — I- -h m+2 h n — I- -h m h n 

7 ftl/i2 dh 2 + 2 hl h2 dh 2 + 2 hlh2 dh% 

m{m - N) n(n - N) f N — I 



hfh% 



— m 



Mi 

E 



ho 



h 2 - Zi 



M 2 

E 

k=l 



ho 



h 2 - w k 



n 



hi 



M 2 , 

+ Ett A 

' hi — 



Mi 

y 

1 hi - Zi f-' hi - w k 

, i=l k=l 



+ 



l-E- 



ft? 



Mi Mi 

EE 



4e ft? 



M 2 



2 (^1 - w fc)(/ii - Wi) 



+ {/ii <-> h 2 } 



(60) 



In this term, the restriction of the allowed momentum eigenvalues (|49|) follows from the 
.B-series as in the one-holon case. 

The other charge kinetic terms are treated by using the fact that the two-holon wave 
function can be expressed by either pairs of color variables, as is shown in App. O The 
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terms involving green variables are rewritten in terms of the z^s and Wk's by the identities 
given in App. |E| Thus we finally deduce for the sum of the three charge kinetic terms 



N 



C aa C t/3 



;<Pr. 



[zf, w k ] hi, h 2 ] 



N 2 - 1 3 / *n 3 , 
1 — m[m — JM) H — n[n — is) 

— 2, 2 



+(n + m)(N-2)- 1C\ - 2C 2 - (m - n) 



N - 1 hx 



Mi , M 2 

ft. _ 7 , 



.4 = 1 fc = l 



hi 



hi 



Mi 

E 



fc 

Mi 



f (/ii - w k ) 2 hx - h 2 V -f-f h - Zi hx - w k 

1 \ 1=1 k=l 



hi + /i 2 
fti - h 2 



M 2 

E 



+ 



hi 



hi 



* (hx-Zi) 2 j-^j (hx- Zi)(hx- Zj) 



M 2 

+E 

k^l 



hi 



+ 2- 



hl 



(hx - w k ) (h x - w{) (hx-h 2 ) 



k=l 

Mi M 2 

i^tri {hx-Zi)(hx-w k ) 
f {hx <-> h 2 } 



ho 
rrim i 



(61) 



where the constants Ci and C 2 are defined as above. 

As can be readily verified, all non-diagonal terms cancel. Summing up the diagonal 
contributions, we obtain the action of f/su(3) on the auxiliary wave functions <p mn , 



ft 



2tt 2 



SU(3) ¥W 



72 (-40 + 33iV - N 3 ) + ^m(m - iV) + - N) 



+ (n + m)(N-2) + 2 



+ ^hx + h 2 



(hx - ^ 



(m — n) 



/i 2 — h 2 



(62) 



<£>r. 



Using (jSHJ), we thus deduce 



#SU(3) *mn 



36 V N N 2 J 



4 ^ 



+ 



3tt2 

'AT2 

2tt 2 



iV+l\ / iV + l\ ;//-// 
m — ) m + ( n — ) n + 



ho 



N 2 



(m — n) ^ 



ho 

m—£,n+£i 



(63) 



=i 



where we have used (x m y n — x n y 

x~y \ & 



x +y („m„,n r "ii'"\ o V" " )•"■'-■' ,,»-!-' ( x m y n + x n y m ) and [J denotes 

the floor function, i.e., [^J is the largest integer / < x. First, note that the action of the 



Hamiltonian on * s trigonal, i.e., the "scattering" in the last line is only to smaller values 
of m — n. Second, shows that the states \I/^ n form a non-orthogonal set, out of which 
we can construct an orthogonal basis of eigenfunctions as it is shown in the following. 
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3. Energy eigenstates 



Using the Ansatz 

ft)= E «rKu*<>, (64) 

for the diagonalization of (|63|). we obtain the recursion relation 



mn 



1 £_1 

T (n-m-2i)ar, 0^ = 1, (65) 



3^ + m-n-|) , () 

which defines the two-holon energy eigenstates (|64|). The corresponding energies are given 



by 



7T / „ a 4 \ 3?r 2 



iV+l\ / N+l\ m-n 
m m + \ n n H 



(66) 



36 V N N 2 ) N 2 [ V 3 

where the momentum quantum numbers are restricted to the interval (j4T?j) and the total 
momentum is given by (|51|). 

The two-holon energies can be rewritten using the one-holon dispersion (|21j) as 

E h Z = Eo-^ 2 + e h °(^°) + 6 ho W°), (67) 

where we have introduced single-holon momenta according to 

ho 2vr 2tt hn 2-k 2-k ( l\ 

We will discuss the physical interpretation of this assignment in Section HVI 



III. SU(n) KURAMOTO YOKOYAMA MODEL 

In this section we extend our investigations to the SU(n) KYM. We will concentrate 
on stating the results and make only short remarks on the calculation, since the decisive 
methods were already discussed in detail for the SU(3) case. 



A. Hamiltonian 



Consider an underdoped chain with at most one particle per lattice site carrying an inter- 
nal SU(n) quantum number which transforms according to the fundamental representation 
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G ■ 



n of SU(n). Starting from the general expression for the SU(n) KYM, the Hamiltonian 
can be rewritten as 

27T 2 1 1 \ - / + 4- \ 1 x - n a na 1 

(69) 

where the summation index cr runs over all flavors 1, . . . ,n, and the Gutzwiller projector 
Pg enforces at most single occupancy on all lattice sites. The model possesses an SU(l|n) 
symmetry generated by the traceless parts of the operators J ab = a) aa a abi where a aa 
annihilates a particle of flavor a at site rj a , as well as a super- Yangian symmetry jis| . 

B. Vacuum state 

We first consider the state containing no excitations. We use a polarized state of particles 
of flavor n as reference state and label the coordinates of the particles of flavor a, 1 < a < 
n — 1, by zf , 1 < i < M a . It can be shown that the states with wave functions |21[ 

n-l M a n—1 Ma- M T n-l M a 

= n n« - 4) 2 n n n« - *d n n < ( 7 °) 

<t=1 i<j <r<r j=l j' = l <t=1 i=l 

constitute exact eigenstates j^] of the Hamiltonian (J69j) . For N = nM, M a = M, i.e., at 
one nth filling, (J70*Jl is the ground state of (jo^jl with energy 

7T 2 fn -2 2ra- 1\ 

e °=--2\— n+ — > < 7i) 

The momentum is p = (n — l)7rM mod 27r, i.e., p = for n odd and p = or p = n 
otherwise. 



C. Spinon excitations 



For iV = nM — 1, localized SU(n) spinons are represented by the wave function 







Mi 

n P [<]=n(^-*9*o[<], (72) 

1=1 

where Mi = M — 1 and M2 = . . . = M n _i = M. The spinons transform according to the 
representation ft under SU(n) transformations. Momentum eigenstates are constructed via 
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Fourier transformation, the spinon momenta are given by 

71 — 1 



P? 



n-1 2vr 

TriV- — 

n N 



v + 



In 



mod 27r, 



(73) 



where the momentum quantum number v is restricted to < v < M\. The momenta (|73|) 
fill the interval [— -, -1 for n even and M odd, or the interval W — —, 7T + -1 otherwise. The 

Ln'nJ ' L n ' n' 

one-spinon energies are given by 

n 2 - 1 vr 2 



= E + 



12n iV 2 



(^ P )- 



(74) 



with 



6 S P( P ) 



n / tt 



A\n 2 

'2 



V 



(p 



7T 



if n even and M odd, 
otherwise. 



(75) 



n n 
, 4 \n 2 

SU(ra) spinons obey fractional statistics, the statistical parameter between spin-polarized 
spinons is given by g = (n — l)/n. 



D. One-Holon excitations 

For iV = nM + 1 one can show by a straight-forward generalization of the SU(3) calcu- 
lations that the one-holon states represented by the wave functions 

n— 1 M a 

vI/JV; h] = h» n II^ - z °) *o[«1, (76) 

(7=1 1=1 

are eigenstates of the Hamiltonian (|69|). In (j76jl . /i denotes the holon coordinate, Mi = . . . = 
M n _! = M, and the momentum quantum number \i is restricted to 

N + n- 1 . x 

< u < — . 77 

n 

The state corresponding to 1)69)1 is constructed in analogy to (jlfij) . The one-holon momenta 
are given by 

hn n — 1 2it / n — 1\ , , . 

# = — ' w + iv("- -s") mod2 "' (78) 

which fill the interval [— -, -l for n even and M odd, or the interval \ir — -, 7r + -1 otherwise 
(either n odd or M even or both). The one-holon energies are 
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(a) n even 

e ho (p) . M odd 







(b) n odd 
e ho { P ) 



M even 



2tt 









2ir V 



FIG. 1: SU(n) holon dispersion, a) n even. The allowed momenta fill the interval [— ~, ^] for M 
odd and \ir — —, tt + — 1 for M even, b) n odd. The allowed momenta fill the interval W — — ,tt + — 1. 



with the single-holon dispersion (see Fig. |TJ) 



6 h °( P ) 





^7T 2 










I 





- (P - 71") 



if n even and M odd, 



otherwise. 



(80) 



E. Two-Holon excitations 

Consider a chain with N = nM + 2 lattice sites. The two-holon momentum eigenstates 
are represented by the wave function 

n— 1 M CT 

/n, h 2 ] = (h - h 2 ){h»h» + h\b%) J] JI^i - z t) (h2 - zf) w k ], (81) 

CT=1 1=1 

where the momentum quantum numbers // and v are restricted to 

N + n - 2 

< i/ < At < . 82 

n 

The total momentum is given by 

, n-1 2tt / n-2\ 

# = — + iv (^ + Z/ ~^J mod27r ' (83) 

As in the SU(3) case, the momentum eigenstates (|81j) form a non-orthogonal basis. The 
two-holon energy eigenstates are obtained using the Ansatz 

A=0 
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where the recursion relation for the coefficients a^T's is found to be 



A-l 



flU 



7iA(A + jz-i/-±) ^ 
The two-holon energies are given by 



IT - ^ - 2k )< 



, u 



1. 



(85) 



ho 



7T 



12n 

r 2 



f(n - 2) A + (2n 2 - 13n + 24)^ - 4(n 2 - 6n + 8)-^ 



+ 



n7r 



AT2 



A + n-2\ / N + n-2' 
H //+ \v I v 



N 2 

\L — V 



, v , (86) 

n J \ n J n 

Using the single-holon dispersions (|80jh the energy eigenvalues of ()86|) can be rewritten as 

where we have introduced single-holon momenta according to 

7T 2-7T / 72 — 3\ . 7T 27T / 71—1 



2n 



7T 27T 



2n 



(88) 



and restricted ourselves to momenta — ~ < < < ^ for simplicity. 



IV. FRACTIONAL STATISTICS 

Fractional statistics in one dimension was originally introduced by Haldane \^ in terms 
of non-trivial state counting rules. Recently, it was realized that the fractional statistics of 



spinons and holons in the KYM mani 



ests itself also in specific quantization rules for the 



individual spinon and holon momenta 17, 18, 23]. Here we apply this interpretation to the 
holon excitations of the SU(n) KYM. 

First, consider holons in the SU(3) KYM. As we have seen in (jfi7|) . the two-holon energies 
are simply given by the sum of the kinetic energies of the individual holons (and the ground 
state energy). This shows that the holons in the SU(3) KYM are free, which is supported 
by conclusions drawn from the asymptotic Bethe Ansatz j2J]. Furthermore, the momentum 
spacing between the individual holon momenta in (}68|) is 

2tt fl 

Vrn ~ Pn ~ ^- 

which reflects the fractional statistics of the holons with statistical parameter g — 1/3. This 
result is consistent with conclusions reached by Kuramoto and Kato 0, 0] from thermody- 
namics, and by Arikawa, Yamamoto, Saiga, and Kuramoto ,25] from the charge dynamics 
of the model. 



+ 



eN , 



(89) 
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For holons in the SU(n) KYM the situation is similar. From (|57jl we deduce that the 
holons are free, whereas the momentum spacings 

Pm-Pn = ^ , ^N , (90) 



N \n 

obtained from (J88|) show that holons in the SU(n) KYM obey fractional statistics with 
statistical parameter g — 1/n. 

Derived in the context of the KYM, this result has implications for SU(n) spin chains 
in general. In the KYM, where the holons are free in the sense that they only interact 
through their fractional statistics, the individual holon momenta are good quantum numbers. 
They assume fractionally spaced values, which for two holons are given by (|9U|). As the 
statistics of the holons is a quantum invariant and as such independent of the details of the 
model, the fractional spacings are of universal validity as well. If we were to supplement 
the KYM by a potential interaction between the holons, this interaction would introduce 
scattering matrix elements between the exact eigenstates we obtained and labeled according 
to their fractionally spaced single particle momenta. These momenta would hence no longer 
constitute good quantum numbers. The new eigenstates would be superpositions of states 
with different single particle momenta, which individually, however, would still possess the 
fractionally shifted values. The effect of the interaction would hence be to turn the integer 
£ on the right-hand side of (J9*U|) into a superposition of integers, while leaving the fractional 
momentum spacing 2n/Nn unchanged. 

Note that regardless of n, the sum of the statistical parameters of spinons and holons 
always equals the fermionic value one, 

n — 1 1 

9sp + 9h = h - = 1, (91) 

n n 

a result consistent with the concept of spin-charge separation characteristic of these models. 

Finally, as models with SU(n) symmetry in general are frequently studied because of 
simplifying features, it is suggestive to ask whether the large-n limit deserves special atten- 
tion in the model we have studied here as well. Briefly, the answer is no. No part of our 
calculation simplifies in this limit, as we obtain terms similar to the ones encountered above 
regardless of the value of n. In the limit n — > oo, g — > implies that the exclusion statistics 
between holons tends towards bosons. This does not mean, however, that the holons in this 
limit behave like free bosons, but rather that their momentum spacings shrink with the nth 
part of the Brillouin zone they are confined to. 
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V. CONCLUSIONS 



In conclusion, we have constructed the explicit wave functions of the one- and two- 
holon excitations of the SU(n) KYM and derived their exact energies. The holons are non- 
interacting or free, but obey fractional statistics with parameter g = 1/n, which manifests 
itself in the quantization of the single holon momenta, which is a general feature of fractional 
charge excitations in SU(n) spin chains. 
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APPENDIX A: GELL-MANN MATRICES 

The Gell-Mann matrices are explicitly given by 
/ n i n \ 



A 1 



1 

1 




A 



A 7 



1 \ 



1 y 


0-i 

i 



A 



^0 -i 0^ 
i 




( o o -%\ 




% 



A a 



A 



\ 



/ 



V 



i o o 

1 
0-2 



^ 1 
0-10 
^0 

^0 ^ 
1 
1 



\ 



They are normalized as tr (A a A 6 ) = 25 a t and satisfy the commutation relations [A a , A b ] = 
2f abc X c . The structure constants f abc are totally antisymmetric and obey Jacobi's identity 



jabc jcde _|_ jbdc jcae _|_ ydac y 



bdc rcae 



cdac rcbe 



0. 



Explicitly, the non-vanishing structure constants are given by f 123 = i, f 1A7 = f 2m = f 257 = 
y-345 = _yi56 = _j367 = / 458 = f 78 = iy/3/2, and 45 others obtained by permutations 
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of the indices. 

The SU(3) spin operators can be expressed in terms of the colorflip operators and the 
charge occupation operator as 

3 



0=1 



APPENDIX B: USEFUL FORMULAS 



For derivations see for example |4l Il3|. 
1. 



N N 

tf = l, £C = tf«W, J{ri a = {-Xf-\ (Bl) 

o=l a=l 



3. 



a=l 



A' 



A 1 - m N + l 



(B2) 



N 



_ 1 M 

T[(Vp - Va) = lim = -• (B4) 



V — !z— = — m, 1 < m < N. (B5) 

^ ^ - 1 2 

APPENDIX C: REPRESENTATION OF WAVE FUNCTIONS 

It is shown that the wave functions can, up to a minus sign, be expressed by any two sets 
of color variables. First, the wave function of the vacuum state ((7J) can be rewritten using 
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green (u) variables as 

Mi M 2 Mi M 2 Mi M 2 

Wofowfc] = (-i) Mi ^n^~^)n^ fc_ ^) 2 nii(' 2i_ ^)n^ii u;fe 

ij^j k<l i=l k=l i=l k=l 



Mi -l 



-I) {-L) h 3 - - , (I I ] 



nsnso*.-* 

where we have used ()B4|) . Accordingly, if we express \l/ i n terms of green and red variables, 
we find 

M 3 M 2 Af 3 M 2 M s M 2 

*o[u,;w k ] = (-l) Ma ^]J(u s -u t )]J(w k -wi) 2 Y[Y[(u s -w k )Ylu s Ylw k 

s^t k<l s=l fc=l s=l fc=l 

{ } { } { } iffiiffid*.-*) 

= (-l) M ^ [ Zi ;w k }, (C2) 



where we again used (|B4|) . and finally set Mi = M2 = M3 = M. 

The same line of argument can be applied to the one-holon wave functions ()17j) 

UfMh - w k ) Uk=i Uti K - m) 2 N^ 

k<l 

nfinSK- 

whereas starting with green and red variables yields 



^[ Zl ;w k ;h] = (-i)*^(-i)«Mh n — „ Ml ~ Ma , kh ; , (C3) 



*%[u a ;w k ;h] 

Ut%(h ~ v, h ) n£i v>k lit K - n) 2 N M * 

= ( _ 1) M 3 ^ ( _ 1) Af 1 Af 3( _ 1) M 2 M 3/l n fc<2 

= (-1)^*51^,; to*; fc]. (C4) 

In the same way we find for the two-holon wave functions (|4*H|) 

fe; *n, fta] = (-l)^*^.; Wfe ; & lf fc]. (C5) 

Thus, the holon wave functions can be expressed by any two sets of color indices. All 
statements generalize to SU(n). 
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APPENDIX D: S-SERIES 



The 5-series is defined as 

N-l 



a=l 

where < £ < 2(N - 1) /3. Now, B™ equals ijEfij). = m - (JV - l)/2 for < m < iV by 
(jEffl). and B? = 1 for < m < N + 1 by (plj) . Furthermore, for 3 < £ < 2(N - l)/3 we 
find 

{N + 2 
0, for < m < , 
f £-2 \ N + 2 ( D2 ) 

Proof: 



a=l fc=0 ^ ' 

= - E( £ : 2 )(- i )"f i -E"« +, ' +i ) 

fc=0 \ ' V a=l / 

fc=0 ^ ' 

Thus, for < m < (iV + 2)/3, I?™ vanishes, as the sums of the binomial coefficients of even 
sites and odd sites equal each other. For (N + 2)/3 < m, however, B™ ^ 0, and thus the 
Taylor expansion appearing in the calculations of the charge kinetic terms contains higher 
order derivatives. 

The remaining constants are deduced from An = B\, C\ = B^" 1 , and C2 = —B^ 1 . 
APPENDIX E: DERIVATIVE IDENTITIES 

If one ho Ion is present, we use for the simplification of the charge kinetic terms 

M 3 , N , Mi , Mi , 

^ — v ft ^ — > fl ^ — v it ^ — \ ll 

' h — fi, / -— ' h — rir, ' /?, — z,- ^— ' h — wi. 

s=l s a^h la i=l 1 k=l K 

AT 1 Ml J, Mz U 

~~ ^ h- Zi~ ^ h-Wk 

i=l 1 k=l K 



\Z-k-2 k 
la 
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as well as 



M 3 h2 M 3 h2 M 3 ^ 



^ — V lb V "\ lb V -V 

j£ (h-u s )(h-ut) ~ ^(h- u a )(h -u t )~j^(h- u s f 



c f^Jh c " 2h f. h* 

M 2 h2 Mi M 2 ^ 

+ ^(h- Wk )(h- Wl ) + ^ 1 ^ l (h-z i )(h- Wk y (E2) 

with the constants C\ and C 2 as above. 

For the two-holon case, we apply the identity 

y hj = _ c , ^ h\ ^ hj h\ 

^ (h - u^h -u t ) 2 + ^ (h - Zt y + ^ (h, - Wk y + (h - h 2 y 

— yj^ — tAt 

^ Aii - ^ f^ht-Wk h- h 2 

i=l k=l 



X 



Mi , M 2 , , 

\ ^ Ail \^ ^1 ^1 



hi- Zj h- wi hi - h 2 

3=1 1=1 

h\ ^ hj 



h 2 

-c 2 + Cj + g { j li _ z .^ hi _ Zj) + _ _ WJ ) 

A 2fr 2 ^ 2fr 2 2fr 2 2fei 

Mi M 2 2 Mi M 2 

+ 2^ 2_^ _ rvr ~r\ - °i 2^ irr ~ " Ul 2^ 



tl £^ C*l - Zi) (^1 ~ ~^ ^1 - ^ ^1 - W k 

| 2fei Ay fei | yv fei \ 

and the similar result for <-> h 2 . All identities presented above directly generalize to 
SU(ra). 



[*] Present address: Rudolf Peierls Centre for Theoretical Physics, University Oxford, Oxford 
0X1 3NP, UK. 
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